Clusters of calcium release channels harness the Ising phase transition to confine their elementary intracellular signals.

Intracellular Ca signals represent a universal mechanism of cell function. Messages carried by Ca are local, rapid, and powerful enough to be delivered over the thermal noise. A higher signal-to-noise ratio is achieved by a cooperative action of Ca release channels such as IP3 receptors or ryanodine receptors arranged in clusters (release units) containing a few to several hundred release channels. The channels synchronize their openings via Ca-induced Ca release, generating high-amplitude local Ca signals known as puffs in neurons and sparks in muscle cells. Despite the positive feedback nature of the activation, Ca signals are strictly confined in time and space by an unexplained termination mechanism. Here we show that the collective transition of release channels from an open to a closed state is identical to the phase transition associated with the reversal of magnetic field in an Ising ferromagnet. Our simple quantitative criterion closely predicts the Ca store depletion level required for spark termination for each cluster size. We further formulate exact requirements that a cluster of release channels should satisfy in any cell type for our mapping to the Ising model and the associated formula to remain valid. Thus, we describe deterministically the behavior of a system on a coarser scale (release unit) that is random on a finer scale (release channels), bridging the gap between scales. Our results provide exact mapping of a nanoscale biological signaling model to an interacting particle system in statistical physics, making the extensive mathematical apparatus available to quantitative biology

Mechanisms of Calcium Leak from Cardiac Sarcoplasmic Reticulum Revealed by Statistical Mechanics

Heart muscle contraction is normally activated by a synchronized Ca release from sarcoplasmic reticulum (SR), a major intracellular Ca store. However, under abnormal conditions Ca leaks from the SR, decreasing heart contraction amplitude and increasing risk of life-threatening arrhythmia. The mechanisms and regimes of SR operation generating the abnormal Ca leak remain unclear. Here we employed both numerical and analytical modeling to get mechanistic insights into the emergent Ca leak phenomenon. Our numerical simulations using a detailed realistic model of Ca release unit (CRU) reveal sharp transitions resulting in Ca leak. The emergence of leak is closely mapped mathematically to the Ising model from statistical mechanics. The system steady-state behavior is determined by two aggregate parameters: the analogues of magnetic field ($h$) and the inverse temperature ($\beta$) in the Ising model, for which we have explicit formulas in terms of SR Ca and release channel opening/closing rates. The classification of leak regimes takes the shape of a phase $\beta$-$h$ diagram, with the regime boundaries occurring at $h$=0 and a critical value of $\beta$ ($\beta*$) which we estimate using a classical Ising model and mean field theory. Our theory predicts that a synchronized Ca leak will occur when $h$>0 and $\beta>\beta*$ and a disordered leak occurs when $\beta<\beta*$ and $h$ is not too negative. The disorder leak is distinguished from synchronized leak (in long-lasting sparks) by larger Peierls contour lengths, an output parameter reflecting degree of disorder. Thus, in addition to our detailed numerical model approach we also offer an instantaneous computational tool using analytical formulas of the Ising model for respective RyR parameters and SR Ca load that describe and classify phase transitions and leak emergence.Comment: 20 pages, 6 figures, supplemental materia

Forehead Skin Blood Flow in Normal Neonates during Active and Quiet Sleep, Measured with a Diode Laser Doppler Instrument

Changes in forehead skin blood flow during active and quiet sleep were determined in 16 healthy neonates using a recently developed semi-conductor laser Doppler flow meter without light conducting fibres. Measurements were carried out at a postnatal age varying from 5 hours to 7 days. The two sleep states could be distinguished in 17 recordings. The mean skin blood flow values during active sleep were significantly higher (p<0.01) than those during quiet sleep, the mean increase being 28.1%. The variability of the flow signal, expressed as the coefficient of variation, changed significantly from 23.1% during active sleep to 18.2% during quiet sleep